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We investigate the influence of degeneracies of the conduction band and the /-orbital on the sta- 
bility of ferromagnetism in the periodic Anderson model. To this end we calculate the temperature 
dependence of the inverse susceptibility for different degeneracies D C ,D^ and conduction electron 
densities n c within the dynamical mean-field theory. A strong increase of the Curie temperature 
T c with the degeneracy of the localized /-level is found. For D c < a simple ansatz based 
on a mean-field treatment of the RKKY interaction is shown to imply a scaling behavior of T c as a 
function of the conduction electron density per band which is well obeyed by the numerical results. 
In particular, T c is found to have a maximum at n c /D c « 0.3. 

PACS numbers: 71.10.Fd, 75.10.-b, 75.20.Hr 
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I. INTRODUCTION 

The periodic Anderson model (PAM) is the mini- 
mal model for the investigation of interacting, local- 
ized electrons (e.g., /-electrons) hybridizing with a 
non-interacting band of conduction electrons as in the 
rare earths and actinides. Indeed, the PAM can ex- 
plain many of the characteristic properties of heavy 
fermion systems^ intermediate valence materials^* 2 - and 
Kondo insulators^ While it is well known that the 
PAM can account for the long-range antiferromagnetic 
correlations at and near half-filling of the conduction 
band, ferromagnetic phases of this model have received 
less attention. Nevertheless, ferromagnetism induced 
by localized^ && £ & £ and even slightly delocalizecU^ /- 
electrons is now known to be a generic property of 
this model away from half-filling. In fact, there ex- 
ist a number of /-electron materials with ferromag- 
netic phases, e.g., UCu 2 Si 2f ±i UCu 2 Ge 2f ii CeRh 3 B 2 r2 
YbNiSn^ URhSijii URhGe,±i and CeRuPO^ Since 
the ions providing the /-electrons in these systems are 
rather far apart, their hybridization with conduction elec- 
trons must be considered essential for the stabilization of 
ferromagnetism. This makes the PAM an appropriate 
model for the investigation of these systems. 

The conventional PAM does not take into account any 
degeneracies of the electrons, e.g., the fact that with- 
out crystal field splitting the /-orbitals have a seven-fold 
degeneracy. This degeneracy can be expected to influ- 
ence the magnetic properties of materials. For example, 
in the case of the Hubbard model band-degeneracy is 
known to greatly enhance ferromagnetism due to pres- 
ence of Hund's rule couplingS i 16 ' 17 i 18 Therefore, a more 
realistic study of ferromagnetism in the PAM should ex- 
plicitly consider band and orbital degeneracies of the 
electrons^ In this paper, we present a detailed study of 
the influence of band degeneracy of the conduction band 
and/or orbital degeneracy of the localized levels on the 



stability of ferromagnetism in the PAM for different con- 
duction electron densities. For this purpose we employ 
the dynamical mean-field theory (DMFT) with quantum 
Monte-Carlo (QMC) as the impurity solver. In partic- 
ular, we show that, and explain why, T c increases with 
the degeneracy of the /-level. The paper is organized as 
follows. In Sec. [IT] we introduce the model and the calcu- 
lation method in detail. The result of the magnetization, 
magnetic susceptibility, and phase diagrams in the PAM 
with degeneracy are presented and discussed in Sec. IIIII 
A summary is presented in Sec. IIV1 



II. MODEL AND METHOD OF CALCULATION 



The multi-orbital PAM investigated here has the form 
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where c i[cr (fj mr7 ) creates conduction (localized) electrons 
with orbital I (m) and spin a at site i. Furthermore, tf is 
the hopping parameter of the conduction electron, while 
ef and denote the center of the conduction band and 
the energy level of the /-electrons, respectively. The or- 
bital index can take the values 1 = 1,..., D c for the con- 
duction band and m — 1, . . . , D* for the /-level. In this 
paper, we consider only degenerate bands, i.e., = t c , 
el = e c , and = £ ^ ■ Furthermore, U and U' are the 
intra- and inter-orbital Coulomb repulsion, respectively, 
F is the Ising component of the Hund's rule coupling, 
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and V represents the hybridization between the conduc- 
tion band and the localized orbital. For the conduc- 
tion band we assume a Bethe density of states (DOS) 
with p c (s) — (l/27r)vi — e 2 per spin. The bandwidth of 
the conduction band (W — 4) defines the energy scale. 
We note that the magnetic properties of the PAM do 
not seem to depend on the specific shape of the non- 
interacting DOS, since a Gaussian DOS gives qualita- 
tively similar results^ In this paper we fix the interac- 
tions, the hybridization, and the chemical potential at 
the values U = 1.5, U' = 1.1, F = 0.2, V = 0.6, and 
fj, = s f + U/2 + (D f - 1)(U' - F/2), respectively. For 
the latter choice of /i, the /-level is approximately half- 
filled^ thus permitting local moments to be formed. 

In the local moment regime of the PAM ferromag- 
netism is due to the magnetic interaction between 
the localized magnetic moments mediated by the non- 
interacting conduction electrons. While a strong hy- 
bridization increases the magnetic interaction, it simul- 
taneously weakens the magnetic moments. For this rea- 
son, static mean-field theories or perturbative approaches 
do not describe ferromagnetic solutions of the PAM ad- 
equately. Here the DMFT 2 — has proved to be a reliable 
non-perturbative investigation scheme. The DMFT is 
based on the limit of infinite dimensions of the lattice, 
where the many-body problem reduces to a local one and 
can, in principle, be solved exactly^ 

The structure of the DMFT self-consistency equations 
for the PAM is the same as that for the Hubbard model, 
except for the calculation of the lattice Green function 
through the k-integrated Dyson equation, 25 which now 
reads 
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Here ~E ma (uj n ) is the self-energy of the /-electron with or- 
bital m and spin a, and Nk is the number of k points in 
the summation. The Green functions are represented as a 
function of the Matsubara frequency u> n = (2n + l)7r//3, 
where f3 — X/ksT. Then, the bath Green function is 
determined by the self-consistency condition G^i^n) = 
S mCT (w„) + [G^ lcr (w„)]~ 1 , and the Green function is ob- 
tained by solving the effective single-impurity problem 
using QMCi2£ This QMC-method involves a time dis- 
cretization At = (3/L and subsequent extrapolation 
At — > 0. In the present work we do not extrapolate in 
every case, and mostly work with At = 0.25, because we 
are mainly interested in the qualitative behavior of T c . 
As will be shown later the value of T c calculated with 
At = 0.25 is within a few percent of the value obtained 
by the extrapolation At — * 0. For small U (smaller than 
the bandwidth) the resulting error is small. As U in- 
creases, the error increases correspondingly. 

The magnetic susceptibility was calculated from the 
two-particle correlation function^ It can also be ob- 
tained from the magnetization in a weak magnetic field. 
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FIG. 1: (a) Inverse susceptibility Xf^ of the PAM as a func- 
tion of temperature for three values of the /-level degeneracy 
— 1 (circles), D ! — 2 (triangles), D* — 3 (squares) and 
two values of the conduction electron density n c = 0.3 (open 
symbols) and n c = 0.5 (filled symbols); the conduction band 
degeneracy is D c = 1 in all cases. Lines were obtained by fit- 
ting to Eq. @. (b) Inverse susceptibility Xf 1 calculated with 
At = 0.25 compared with the extrapolated result (At — > 0) 
for D f = 3 and n c = 0.4. 



We found that both methods give the same results, al- 
though for the same computational time the statistical 
error in the latter method is much larger. 



III. RESULTS AND DISCUSSION 

The numerical results obtained for the inverse ferro- 
magnetic susceptibility Xf 1 of the PAM for values of 
the /-level degeneracy D* — 1,2,3 and electron densi- 
ties n c = 0.3,0.5 are shown in Fig. [Ha). In the non- 
degenerate case [D* — 1), the inverse susceptibility is 
found to be almost perfectly linear for all temperatures. 
A linear behavior of the inverse susceptibility near T c , 
i.e., a Curie- Weiss mean- field behavior, does not come 
unexpected in DMFT, although its validity was so far 
proved only for the single-band Hubbard model. 28 The 
numerical finding that this linear behavior persists up 
to high temperatures, i.e., that the Curie- Weiss behavior 
merges into a Curie law, is also not surprising in view of 
the low values of T c and a saturated magnetization in the 
magnetic phase. 

For increasing degeneracy D* , the inverse susceptibil- 
ity becomes more and more non-linear. As can be seen 
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FIG. 2: Magnetization and inverse susceptibility Xp 1 extrap- 
olated to At — > as a function of temperature for D c — 1, 
D = 3, and n c = 0.4. The dashed line marks a root- 
like disappearance of the magnetization fitted to the last 
two data points in the magnetic phase. The inset shows 
the At dependence of the magnetization for temperatures 
T = 1/26, 1/24, 1/23, and 1/22. 



from Fig. [Tfb) , this is not due to the finite value of At 
since an extrapolation At — > does not change the re- 
sult. A non- linear behavior of the susceptibility is also 
seen in the orbitally degenerate Hubbard models This 
non-linear behavior, i.e., the change of slope with tem- 
perature, is due to the decrease of the effective magnetic 

T c ) with temperature. This 

1). 



moment ^ cff = \J X F {T ■ 

effect is negligible in the non-degenerate case (D* 
By contrast, in the degenerate case the spins of the lo- 
calized electrons are aligned parallel due to the Hund's 
rule coupling near T c (which is here found to be low) 
such that the effective magnetic moment becomes large. 
On the other hand, at high temperatures the effective 
magnetic moment decreases since the electrons become 
independent and the spin of the localized electron takes 
a random orientation on every site, owing to the ther- 
mal fluctuation. Accordingly, the decrease of the effec- 
tive magnetic moment with increasing temperature leads 
to an increase of the slope of the inverse susceptibility. 
This explains the upturn of the inverse susceptibility for 
increasing temperature. 

The non-linearity of x^ 1 has to be taken into account 
in the calculation of T c from the inverse susceptibility. 
Namely, a purely linear fitting of x^ 1 would lead to a 
serious error in the value of T c . We found the inverse 
susceptibility to be well fitted by 
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1 = A(T - T c ) + B(T - T c 



(4) 



The linear behavior is restored again at high tempera- 
tures, i.e., above the Fermi temperature (Curie law). The 
two linear regimes at high temperatures and near T c are 
smoothly connected. 

To check the validity of Eq. ^ and the T c values ob- 
tained from that expression, we calculated the magneti- 
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FIG. 3: Curie temperature as a function of conduction elec- 
tron density n c for different values of the /-level degeneracy 
and conduction band degeneracy D c . 



zation = (n|) — (n|) of the ferromagnetic phase of 
the degenerate PAM. To reduce the error from thermal 
fluctuations in the QMC calculation we performed an 
extrapolation 1/Nmc — > 0, where A?mc is the number of 
Monte-Carlo steps in the QMC routine, to calculate the 
magnetization at each iteration. Then we performed an 
extrapolation At — > 0. As shown in the inset of Fig. [21 
the magnetization depends linearly on At and a spon- 
taneous magnetization begins to appear at temperature 
between T — 1/23 and 1/22. The magnetization and in- 
verse magnetic susceptibility obtained after the extrap- 
olation At — > is presented in Fig. [51 The value of T c 
obtained by fitting the susceptibility to Eq. (J4| is clearly 
consistent with the values inferred from the magnetiza- 
tion (inset of Fig. \2§ and from the root-fit to the magne- 
tization (dashed curve in Fig. [5]) . 

The dependence of the Curie temperature T c on the 
conduction electron density n c is shown in Fig. [5] for 
several values of the /-level (D*) and conduction band 
(D c ) degeneracy. The most remarkable feature is the pro- 
nounced increase of T c with D f . This behavior can be 
explained by the increase of the local magnetic moment. 
Namely, the Hund's rule coupling aligns the local mag- 
netic moments on each site whereby the magnitude of 
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FIG. 4: Rescaled plot of the Curie temperature T c vs. con- 
duction electron density n c for several values of the /-level 
degeneracy D' and conduction band degeneracy D c , with 
D c <D f . 



tron density per band (n° s ). Therefore, when we plot 
T c /{(Df /2)(D f /2 + 1)/[U+ (D f - l)F} 2 } as a function 
of n c /D c , all curves should fall onto the same curve. The 
result is shown in Fig. [3J The deduced scaling behavior 
is seen to be obeyed surprisingly well, given the simplic- 
ity of the assumptions made in the derivation and the 
fact, that the parameter values employed in our inves- 
tigation are not limited to the RKKY regime. Indeed 
T c increases rapidly with until it reaches a maxi- 
mum near n£ ff = 0.3 and then decreases slowly. The 
maximal value of T c increases with Df and is about 
0.05{(Df /2)(Df /2 + l)/[U + {Df - 1)F} 2 }. In accor- 
dance with the mean-field character of the above ansatz 
the curves for D c > 1 are found to fulfill the scaling 
behavior much better than for D c = 1. For D c > £)•', 
which are not plotted in Fig.[¥J the above ansatz does not 
apply. The additional conduction bands seem to impede 
the magnetic ordering rather than mediate it. 



the local moment, S, increases with the degeneracy Df; 
T c then increases accordingly. By contrast, the degen- 
eracy of the conduction band has a different effect: For 
increasing D c the conduction electron density giving the 
maximal T c (i.e., the "optimal" density) also increases 
whereby the ferromagnetic phase expands to higher n c 
values. At the same time T c is hardly affected. For 
D c > D> , however, T c almost always decreases with in- 
creasing D c . 

In order to explain the general dependence of T c on 
the degeneracy shown in Fig. [3] we modify the ansatz 
that we recently introduced in the case of the disor- 
dered PAM without degeneracy^ which has the form 
T c (p) = T°Ff(n - e f )F c (n - e c ). Here the two func- 
tions F' and F c describe the formation of the local /- 
electron moments and the mediation of the magnetic or- 
dering by c-electrons, respectively. Since \x — is fixed, 
(fx — e?) depends only on the degeneracy. Accord- 
ing to the mean-field theory for the RKKY interaction 
in which this ansatz is indeed justified the local mo- 
ment part (F*) is proportional to J 2 S S(S + 1), where 
J ff = — 8V 2 /U e s^- The effective Coulomb interaction 
U cS is defined as U eS = E(n f + l) + E(n f Q - 1) - 2E(n f ), 
where E(uq) is the energy with the /-electron den- 
sity UqM^ For a half-filled /-level, — D f and 
U eS = U + (Df - 1)F. Thus, assuming S = D / /2, 
F? is proportional to (Df /2)(Df /2 + 1)/[U + (Df - 
1)F] 2 . As in the conduction band part, F c (/i — e c ) can 



be expressed by F c (n^ s ) with 



"eff 



n c /D c , because 



fi — e c is uniquely determined by the conduction clcc- 



IV. CONCLUSION 

We computed the ferromagnetic susceptibility of the 
periodic Anderson model (PAM), and from that the 
Curie temperature T c , for different values of the conduc- 
tion band and /-level degeneracies as well as electron 
densities n c within DMFT. The magnetic susceptibility 
was found to deviate from the Curie- Weiss behavior away 
from T c , an effect which becomes stronger with increasing 
degeneracy of the /-level. Our results show that the de- 
generacy affects the ferromagnetic stability of the PAM 
strongly. In particular, for D c < Df the value of T c (i) 
increases with the /-level degeneracy, (ii) approximately 
obeys a scaling law as a function of the conduction elec- 
tron density per band n c /D c , which can be understood 
within a mean-field picture of the RKKY interaction, and 
(iii) becomes maximal for n c /D c w 0.3. This shows that 
realistic investigations of correlated electron materials 
with localized /-levels must employ appropriate gener- 
alizations of the PAM where the multi-orbital character 
of the electrons is taken into account. 
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